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Abstract 

We justify a recently proposed prescription for performing Green Function 
Monte Carlo calculations on systems of lattice fermions, by which one is able 
to avoid the sign problem. We generalize the prescription such that it can 
also be used for problems with hopping terms of different signs. We prove 
that the effective Hamiltonian, used in this method, leads to an upper bound 
for the ground-state energy of the real Hamiltonian, and we illustrate the 
effectiveness of the method on small systems. 
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I. MOTIVATION 



As is well known, exact Monte Carlo methods cannot be applied straightforwardly to 
fermionic systems. In such systems, the sign problem causes great difficulties in obtaining 
sufficient statistical accuracy,0 particularly as the number of quantum particles increases. 
The reason is that, when sampling physical properties in configuration space, one collects 
large positive and negative contributions, due to the fact that a fermion wave function is 
of different sign in different regions of the configuration space. These contributions tend 
to cancel, giving a result that may be exponentially smaller than the positive and negative 
contributions separately. 

Recently, some of us described a method to perform Green Function Monte Carlo 
(GFMC) on a system of fermions on a lattice^, which is an extension of the fixed-node Monte 
Carlo method for continuum problems, developed by Ceperley and Alder!. In this method 
one avoids the sign problem, replacing the original Hamiltonian by an effective Hamiltonian, 
such that one obtains contributions of one sign only in the sampling procedure. The price 
one has to pay is in the fact that the ground-state energy E e g of the effective Hamiltonian 
is in general not the same as the ground-state energy Eq of the original Hamiltonian. It was 
claimed, however, that E e R is a true upper bound for Eq, making the method variational. 

The proof for this upper bound is less obvious than was suggested in Ref. ^ because an 
assumption was used about the form of the ground state of the effective Hamiltonian, which 
is not generally true [see the discussion following Eq. (|PTD]. It is possible, however, to give 
a general proof for the upper bound. In the process of deriving this proof, we found that 
our prescription can be generalized, such that also problems with a Hamiltonian containing 
hopping terms of different signs can be treated by this method. The aim of this paper is 
to give a general proof, illustrate the method on small systems, for which we diagonalize 
both the original and the effective Hamiltonians exactly, and discuss the applicability of the 
method. We do not actually perform Monte Carlo simulations here. 



II. EFFECTIVE HAMILTONIAN 

We work in a configuration space {R}, where each R denotes a configuration of numbered 
fermions on a lattice. In this configuration space, the Hamiltonian 7i of our problem can 
be represented by a real symmetric matrix H with elements (R\H\R'). One is generally 
interested in finding the ground-state energy of this Hamiltonian subject to some symmetry 
constraints, for example that the wavefunction be antisymmetric. We suppose that the 
ground state \ipo) of Ti is reasonably well approximated by a trial state \iPt), which is 
defined through its wave function in all possible configurations: ip^{R) = We 
restrict ourselves to real trial wave functions, because the ground-state wave function can 
be taken real in this problem, and the sign of the trial wave function is one of the key 
ingredients for our method. Complex Hamiltonians and trial functions can be treated with 
the so-called fixed-phase method^. Typical examples of the Hamiltonians considered here are 
the Hubbard Hamiltonian or the Kondo lattice model, and the typical trial wave function 
is a determinant obtained by a mean-field approximation. 
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In general, a trial wave function divides the configuration space into nodal regions. A 
nodal region is a set of configurations in which the trial wave function has the same sign, 
and which are connected via the Hamiltonian. For an antisymmetric wave function, there 
is equivalence between the regions of positive and the regions of negative sign. 

In GFMC with importance sampling, random walkers diffuse and branch through the 
configuration space in a stochastic way, guided by a trial wavefunction. The Hamiltonian 
is used to project out the lowest energy state. The process for a lattice model and its 
mathematical justification is described in more detail in Appendix A and Ref. || In the 
fixed-node approach, one ensures that the contribution of a specific walker is always positive, 
otherwise the negative signs will eventually interfere destructively. For completeness, and 
to indicate the connection with the sign problem in quantum Monte Carlo simulations, we 
expand on this point in Appendix A. If the off-diagonal terms in the Hamiltonian are all 
negative (as in the Hubbard model), a sign change only occurs when a walker goes from one 
nodal region to the other. More generally, a walker could collect an unwanted minus sign if 
there exists a pair of configurations R and R' such that: 

{R\H\R')i; T (R)ij T (R') > 0. (1) 

In order to prevent this from happening, we make an effective Hamiltonian which does not 
have such matrix elements. 

The fixed-node method was developed for the case in which the electron coordinates are 
continuous variables!'!. There, one has to deal with kinetic terms of negative sign only, and 
the nodal surface of a trial wave function is uniquely defined as the set of configurations 
where it vanishes. The fixed-node constraint can be implemented by imposing the boundary 
condition that tj; must vanish on the nodal surface of ifa. In the limit of sufficiently small 
step sizes, we can make sure that Eq. (P is never violated since R and R' become closer 
together and t/'t will vanish. In this way one obtains the lowest energy under the condition 
that the wave function has the same nodal structure as the trial wave function. This energy 
yields an upper bound to the true ground-state energy; in practice, very accurate estimates 
for the ground-state energy of continuum problems can be obtained. 

On a lattice, one has to deal with discrete steps, and one has to treat the hops that 
cause a change of sign in a different way. In our implementation, we replace those unwanted 
hopping terms in the Hamiltonian by diagonal terms, that depend on the ratio of the trial 
wave function in the configurations R and R' . We thus construct an effective Hamiltonian 
H c fi as follows: 

(R\H eS \R'} = (R\H\R f ) (if (R\H\R'}ifj T (R)ifj T (R') < 0) 

= (otherwise) (2) 

are the off-diagonal terms, and the diagonal terms are given by 

(R\H cS \R) = (R\H\R) + (R\V S{ \R). (3) 

The last term in (|3]) is the sign-flip potential at R, which corrects for the contributions of 
the steps left out in H eS . V s f has only diagonal elements, which are defined by 
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(R\V si \R)=Y,(m\R')^^l- (4) 



R' 



MR) 



Here the summation is over all [neighboring] configurations R' of R for which (|]) holds. Note 
that this is a significant extension of the prescription presented in Ref. |2|, where we only 
considered the case that all hopping terms are of negative sign, such that the sign-flipping 
hops would coincide exactly with the hops to a different nodal region. In the general case, 
we prefer to speak about sign-flipping steps instead of nodal-boundary steps, as the latter 
term may cause confusion. 

Clearly, by this prescription, a hop that would induce a sign change is replaced by a 
positive diagonal potential. If instead one used only the truncated Hamiltonian as given 
by (0), with the original diagonal matrix elements (R\H B g\R) = (R\H\R), then the value of 
the wave function at the node would be too high and its energy too low. This was found 
in an earlier attempt to perform fixed-node Monte Carlo on lattice fermions by An and van 
Leeuwen.l 

A somewhat similar procedure, called "model-locality", has been used by Mitas et a/.i 
in continuum problems with a non-local potential that arises from replacing atomic cores 
with pseudopotentials. As in a lattice system, they can not solve the problem of crossing a 
node by making the step size of the walkers continuously smaller, because of the non-local 
potential that connects configurations at finite distances. In their approach the unwanted 
off-diagonal terms are truncated, and replaced by diagonal contributions as in Eq. (f|), but 
with the sum over all R', not just over sign-flip configurations. With the model-locality 
procedure, one does not obtain an upper bound for the ground-state energy. 



III. UPPER BOUND 

We want to show that the prescription given above for TC e s leads to an upper bound for 
the ground-state energy of Ti. In order to do so, we define a truncated Hamiltonian H, tT , 
and a sign-flip Hamiltonian 7i s f, by 

H = Htr + H-si, (5) 

n cS = n tv + v rfj (6) 

where the diagonal elements of H tr are 

(R\H tI \R) = (R\H\R), (7) 

and its off-diagonal elements are given by 

(R\H tT \&) = (R\H eS \R'). (8) 

V s f is the sign-flip potential, for which the matrix elements are given by (f|), and H s { contains 
only the off-diagonal elements of H which are put to zero in the effective Hamiltonian. We 

now take any state 

IVHEWW (9) 

R 
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and we compare its energy with respect to Ti and to Tt e $'- 



AE={iP\(H eS -H)\ij) 
AE can be written explicitly in terms of the matrix elements of V s { and H s f. 



R 



(R\V si \R)ifj(R) - J2(R\H si \R')ip(R') 

R' 



(10) 



(11) 



We rewrite this expression in terms of the matrix elements of H: 



AE = j2mr 

R 



sf 



sf 



T,(R\H\R')^lHR) - T,(R\ H \R'MR') 



R> 



MR) 



R' 



(12) 



In this double summation each pair of configurations R and R' occurs twice. We combine 
these terms and rewrite (O) summation over pairs: 



sf 



AE 



E (R\h\#) 

(R,R') 



MR)f 



MR') 
MR) 



+ 



MR) 
MR') 



- iPiRY^R') - ip(R')*^(R) 



(13) 



Denoting by s(R, R') the sign of the matrix element (R\H\R'), and using the fact that for 
all terms in this summation the condition ([]]) is satisfied, we can finally write AE as 



sf 



AE 



E \(R\h\R') 

(R,R') 



HR) 



\ 



MR') 



MR) 



s(R,R')ip(R') 



\ 



MR) 



MR') 



(14) 



Note that we do not have to worry about configurations R where ip^{R) = 0: they do not 
occur in this summation. Obviously, AE is positive for any wave function ip. Thus the 
ground-state energy of 7^ e ff is an upper bound for the ground-state energy of the original 
Hamiltonian Ti. 

Now the GFMC method can calculate the exact ground-state energy E e s and wavefunc- 
tion ^eff of H c s, without any sign problem. Assuming the trial function ip^ has the correct 
symmetry [for example is antisymmetric], then ip e Q will carry the same symmetry and hence: 
E e s > {il>ett\H\ip e ff) > Eq, where the second inequality follows from the usual variational prin- 
ciple. Hence the fixed-node energy is an upper bound to the true ground-state energy. One 
can easily verify that H.\i/jt) = TCes\ipT)-, an d thus one can be sure that the GFMC procedure 
improves on the energy of the trial wave function: E eS < (ilj T \H eS \ip T ) = (if) T \H\ip T ). 



IV. VARIATION OF THE TRIAL STATE 

Let us consider the situation where we use the exact ground state |^o) of H, with energy 
E , as trial state. Obviously, for the method to be useful, it is desirable that in that case 
the effective Hamiltonian has the same ground-state energy E , and the same ground state 
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|?/>o), as that would make it possible to find the true ground state by varying the trial wave 



function in some way. In Eq. (14) we substitute ipo for ifa- In order to have AE equal to 
zero, each individual term in the summation fll~4] ) has to vanish, thus leading to 



MR') 



MR) 



- s(R, R')ip(R') 



\ 



MR) 



or, 



s(R, R!) 



MR) 



MR') 



MR') 



MR) 
MR') 



(15) 



(16) 



for all sign-flipping pairs (R, R'). This condition is trivially fulfilled for ifi = i/jq. Thus, the 
true ground-state energy can be reached by variation of the trial wave function. One can 
further extend this result to show that as M ~^ Vto the error in the fixed-node energy will 
be second order in the difference, M ~ V'oj with the coefficient positive. 

The original contention in Ref. |2] was that the wave function, obtained through this 
effective Hamiltonian, would have exactly the same ratio at each sign-flipping pair R and 
R' as the trial wave function, i.e. 



M{R) = MR) 
i/MR>) MR') 



(17) 



This, however, is in general not the case, as one can see from considerations about the 
symmetry of the wave function. An example for a small system, which illustrates this point, 
is given in Appendix B. Note that our proof for the upper bound does not rely on the 
assumption (fl7|) , and that the conclusion we put forward in Ref. § about the variational 
principle, remains unchanged. In fact, because the ground state of H e s is found in a much 
less restricted space of states than those satisfying (|17|), the resulting estimate for the ground- 
state energy is much better than was anticipated. 

There is an important difference between the lattice and continuum fixed-node method. 
In the continuum method, it is only the sign of the trial function that matters. If the nodes 
are correctly placed, one will obtain the exact energy regardless of the magnitude of the 
trial function. Clearly this does not hold with the lattice fixed-node procedure: the sign of 
the trial function and the relative magnitudes of the trial function in configurations that are 
connected by a sign flip must be correct. For example, in the continuum the exact result 
would be obtained for a one-dimensional problem since the nodal surface are the coincident 
hyperplanes. One does not nessecarily get the exact result for a 1-d lattice model as one of 
the following examples shows. 



V. ILLUSTRATIONS 

We illustrate the effect of the effective Hamiltonian for the single-band Hubbard model 
by exact calculations on two small systems: a loop of 4 lattice sites on the corners of a 
square, and a graph consisting of 8 points on the corners of a cube. We use the well-known 
Hamiltonian 
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W = -*E <i*c s * + U y En tt n il , (18) 

containing nearest-neighbor hopping with strength t and a local interaction between electrons 
of opposite spin with strength U. We consider the square with 2 electrons with spin up 
and 2 with spin down, and the cube with 4 up and 1 down and with 4 up and 4 down, 
respectively. In all cases, we use the results of self-consistent mean-field calculations as 
trial wave functions, and we compare the mean- field energy (MF), the lowest energy of the 
effective Hamiltonian (FN, for fixed node), and the exact ground-state energy, for different 
values of the interaction parameter U. We use different restrictions on the average number of 
electrons with spin up and down per site, in order to obtain different types of self-consistent 
mean-field wave functions. Writing (n ia ) = (n a ) + (—l) a q ia for the average number of spin-cr 
particles on site i, we denote = by H (homogeneous), and qi a = (— lYq a , with q a a 
constant, by AF (antiferromagnetic, or Neel order favoured). In the case of the cube with 4 
up and 1 down spins [i.e. off half filling], the self-consistent mean-field solution with lowest 
energy turns out to have a symmetry different from both H and AF. As the exact ground 
state of TL as well as of 7^ e fr is not degenerate in these cases, it cannot have broken symmetry 
for (n ia ). Note that, for U = 0, the mean-field approximation yields the exact ground state, 
and we checked that also the fixed-node result equals the exact ground-state energy in that 
case. The results are presented in Table [|J 

As one can see, the fixed-node approach on these small systems yields a significant 
improvement on the upper bound for the ground-state energy, compared to the mean-field 
approximations. One may note the fact that the mean-field wave function with lowest energy 
does in general not give the best fixed-node result. In a real problem, one would want to find 
the best possible trial wave function as input for the fixed-node procedure, and it is clear 
from these results that 'best' does not mean 'having the lowest variational energy' here. The 
sign of the trial wave function and its behavior at the nodal boundary determine how good 
the fixed-node energy will be. 

VI. CONCLUSIONS AND OUTLOOK 

This method can be applied to any lattice model provided that only the Hamiltonian 
and a trial wave function with the proper symmetry are given. It is super-variational, in the 
sense that it always yields an upper bound for the energy, which is the lowest possible value 
consistent with the imposed constraint. By varying the the effective Hamiltonian through 
the trial wave function, in principle the exact ground-state energy can be obtained. 

Note that we have not used the symmetries of the trial state as input for our method. 
This means that it is possible to use this method for models of frustrated spins on a lattice 
and, via the appropriate mappings, for systems of bosons as well, or for excited states which 
are ground states of a given symmetry. In a forthcoming publication, possibilities to do so 
will be presented and discussed. 

Note further that the nodal relaxation method, as described in Ref. ^| for continuum 
problems, is also applicable on the lattice. In this method, one uses the fixed-node approach 
to improve on the trial wave function. When this has been done, one removes the sign-flip 
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constraint, and allows the walkers to move through the whole configuration space. If the 
fixed-node result is close enough to the ground state, one can sample the exact ground-state 
energy before the sign problem destroys the accuracy. 

In the near future, we plan to use this method for Monte Carlo studies on some of the 
systems mentioned above, in order to find more comparisons of the fixed-node approach with 
known results, to check the effectiveness of the method, and to tackle some new problems 
as well. 

The method appears to be also useful for continuum problems, where one has a non-local 
potential. For example one can modify the model-locality approach of Ref. |8] so that it does 
yield an upper bound. Essentially one must allow non-local moves which do not change the 
sign and add terms to the effective Hamiltonian corresponding to discarded moves. 

We finally note that another promising avenue for further development of the conceptual 
basis of our approach is given by the observation by Martin that the use of an effective 
Hamiltonian can be couched in the language of Density Functional Theory.0 This makes it 
possible to apply a number of well-known results for the behavior of the energy functional 
under variation of both the effective Hamiltonian and the trial state |^t}- 



The authors want to thank P.J.H. Denteneer for stimulating discussions. This research 
was supported by the Stichting Fundamenteel Onderzoek der Materie (FOM), which is finan- 
cially supported by the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO). 
DMC was supported by the Institute for Theoretical Physics at the University of California 
at Santa Barbara. 



In this Appendix, we clarify the origin of the sign problem for a specific way of per- 
forming Green Function Monte Carlo on lattice fermions, and we explain how one is able 
to circumvent this problem, using the fixed-node approach. More details of this version of 
GFMC as applied to lattices are given in Ref. |5]. 

In a GFMC simulation one tries to obtain information about the properties of the ground 
state of a given Hamiltonian 7i. Starting from a trial state, one can obtain [a stochastic rep- 
resentation of] the ground state by repeatedly applying a projection [or diffusion] operator. 
On a lattice it is simplest to use an operator that is linear in 7i, and that can be viewed as 
the first-order expansion of an exponential diffusion operator in imaginary time: 



where w is a parameter that should be chosen close to the ground-state energy in order to 
keep the wavefunction normalized. The parameter r is taken small enough to ensure that 
the diagonal terms of this operator are positive. The off-diagonal elements in the matrix 
representation for T are, up to a factor — r, the same as those for Ti. The n-th approximation 
of the ground state is given by 
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APPENDIX A: THE SIGN PROBLEM IN MONTE CARLO 




(Al) 
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I^H^IVt). (A2) 

One can check that, if the trial state has some overlap with the ground state, \ip n ) will 
converge exponentially fast to the ground state for large n. 
The ground-state energy can be calculated as: 

W-tIW) , 
E " - IMPT (A3) 

We rewrite this expression as a summation over paths in configuration space: 
E = EiiE(R n )(^\R n )U?= 1 (Ri\F\Ri- 1 )(I^\'ih) 

En{^T\Rn)Utl(Ri\F\Ri-l)(^T) ' 

where E(R) = (i/)t\H\R) (-0t|-R) _1 is the local energy at R, and 1Z = {R , Ri, R 2 , R n } 
denotes a path in configuration space. 

In a GFMC procedure this expression is sampled stochastically by constructing paths 
TZ in configuration space, and calculating the energy from the contributions of those paths. 
Importance sampling is used to reduce the fluctuation of those paths by modifying T . The 
sign problem arises from the fact since the fermion trial wave function is antisymmetric, its 
sign will vary. Also, the matrix elements (R\F\R') between different R and R' need not be 
always positive. Thus, when performing a random walk to obtain a path TZ, starting from a 
configuration i?o where the trial wave function ^t(-Ro) is of specific sign, one may end up in 
a configuration Rn where the trial function is of the opposite sign, or one may have collected 
an odd number of negative (R\F\R') in the path. For large n, one obtains about as many 
positive as negative contributions; the difference is used to determine the energy. One can 
easily show that the "signal-to-noise" ratio must decrease exponentially in n once negative 
contributions are allowed. Intuitively, it is easy to understand that this will give rise to an 
inaccurate result. In practice, this severely limits the applicability of Quantum Monte Carlo 
methods to fermion problems. 

In the fixed-node approach, one wants to avoid that contributions of different sign can 
be obtained. In order to ensure this, one demands that at every individual step along a 
path, only positive contributions are allowed. Thus, all steps satisfying Eq. ([!]) have to be 
discarded. The prescription for the effective Hamiltonian takes care of this constraint. 

Finally, we remark that this prescription fits very well with the way we perform importance 
sampling. When using the trial wave function as a guiding function for the random walks, 
at any point in the walk one needs to know the value of the trial wave function, and one 
can use this value at the same time for guiding the walks and for the implementation of the 
fixed-node effective Hamiltonian. Note that the summation needed to define the effective 
potential in Eq. (|]) only grows linearly with the size of the system for a Hamiltonian such 
as the Hubbard model. Thus it does not appreciably slow the calculation. 



APPENDIX B: EXAMPLE OF FIXED-NODE PROCEDURE 

In this Appendix, we give an illustration of how the effective Hamiltonian is created, 
and what its effect is, on a very simple small system. All steps can be straightforwardly 
generalized to more complicated systems. 
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Consider the Hamiltonian 

n = -t ]T dc ja (bi) 

on a loop of 4 sites with 2 spinless fermions. We define configurations of labeled fermions 
[M2]) where particle j (1 < j < 2) sits on site ij (1 < ij < 4). We number the sites, as 
follows: 

:. 2 

4 3 • 



A valid [i.e. antisymmetric] fermion wave function if) must satisfy if) ([ij] J = 

The configuration space of this system consists of 12 configurations, and can be depicted as 

follows: 




The lines [or bonds] represent valid hops in this space. The matrix elements ([ij]|if|[M]) of 
the Hamiltonian for this system are —t if there is a bond between [ij] and [kl] [in that case 
% = k or j = I must hold], or otherwise. The ground state of this Hamiltonian is symmetric 
under exchange of the particles, and we have to restrict the wave function explicitly to be 
antisymmetric in order to find a valid fermion wave function. To obtain a Hamiltonian H 
which describes the fermion problem only, we define antisymmetric states [ij], which are 
antisymmetrized combinations of the configurations: 

m = ^ (M - m ■ (B2) 

In this way each pair of configurations [ij] and [ji] produces two states, [ij] and [ji], which 
only differ by their sign. One has the freedom to choose one of these states to obtain only 
one state per pair of configurations, and one can calculate the resulting Hamiltonian for the 

z rii n 2 

= ag(U){[ij]\H\U[kl]), (B3) 

where Hi and H 2 denote permutations of the two particles, sg gives the sign of a permutation, 
and H[kl] is the permutation of [kl] that can be reached by one hop from [ij], such that 
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([r/']|i/|n[A;/]) = —t. We can again denote the Hamiltonian in a picture, representing matrix 
elements —t by thin lines and +t by thick lines [we choose the [ij] with i < j; other choices 
give different pictures but the same results]: 



[12] [14] [23] [34] 



This structure fully contains the antisymmetry, and the corresponding Hamiltonian gives 
all information there is on the fermion problem. Its ground-state is degenerate, with energy 
— 2t, and possible ground states are 



and 



+ 



+ 



+ 



(B4) 



Wo) 



|[12] + [13] - [24] 



(B5) 



It is easy to generalize this procedure for any system of lattice fermions. 

Let us now consider a trial state, and calculate the effective Hamiltonian according to 
our fixed-node prescription. The trial wave function defines the nodal regions through its 
sign in all states, and, because we are working with negative hopping terms, the sign-flip 
constraint reduces to sign changes of the wave function only. We take a very simple trial 
state: 



i 



j= | [12] + [13] + [14] + [23] + [24] + [34 



(B6) 



purposely chosen such that we only have to slightly adapt the previous picture to denote 
the effective Hamiltonian: 




Here the thin lines are still matrix elements —t. The thick lines have been cut [we do not 
allow these hops in the effective Hamiltonian] and replaced by arrows, indicating diagonal 
matrix elements, which in this simple case all become +t, because we have chosen equal 
weights for all the states in the trial wave function. The [nondegenerate] ground state of 
this effective Hamiltonian is 

\^ s ) = |0.165([12] + [34]) + 

+ 0.448Q13] +24]) + 0.523([14] + [23])) (B7) 

with energy -1.709t. Note that, e.g., the states [12] and [24] do not have the same wave 
function in this ground state, while they do in the trial state. As one could have expected 
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from symmetry considerations, the wave function is the same in states that have an equiv- 
alent position in the picture, i.e. occur symmetrically in the effective Hamiltonian. States 
that are connected via the boundary do not in general have such symmetry, and thus there 
is no reason to expect that they would obey (|17D . Note also that the energy of the effective 
ground state is above the ground-state energy of the true problem, as it should be according 
to our proof that it is an upper bound for that energy. 
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6 




-6.0701 


-7.2424 


-7.2533 




10 


TT 

H 


-4 


-6.8400 


-6.8442 




10 


A TT 1 

At 


A on ^ 1 
-4.Z001 


-6. (4/0 


-O.o44z 




10 




-5.3271 


-6.7637 


-6.8442 


cube 







-12 


-12 


-12 




1 


H 


-10 


-10.1148 


-10.1188 




2.5 


H 


-7 


-7.7257 


-7.7510 




2.5 


AF 


-7.0061 


-7.6942 


-7.7510 




10 


H 


8 


-2.6597 


-2.8652 




10 


AF 


-2.3113 


-2.6382 


-2.8652 



TABLE I. Comparison of the energies obtained for three different systems by means of 
self-consistent mean-field (MF), fixed-node (FN), and exact calculations^. All values are given 
in units of t. 
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